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We describe the zero-temperature phase diagram of a model of bosons, occupying sites of a linear 
chain, which obey a hard-exclusion constraint: any two nearest-neighbor sites may have at most 
one boson. A special case of our model was recently proposed as a description of a "tilted" Mott 
insulator of atoms trapped in an optical lattice. Our quantum Hamiltonian is shown to generate the 
transfer matrix of Baxter's hard-square model. Aided by exact solutions of a number of special cases, 
and by numerical studies, we obtain a phase diagram containing states with long-range density-wave 
order with period 2 and period 3, and also a floating incommensurate phase. Critical theories for 
the various quantum phase transitions are presented. As a byproduct, we show how to compute the 
Luttinger parameter in integrable theories with hard-exclusion constraints. 



I. INTRODUCTION 

In recent years, the study of quantum models with 
multiple competing ground states has emerged as an im- 
portant theme in the study of strongly-correlated many- 
body quantum systems. For example, in the cuprates it 
is clear that states with density-wave order at a variety 
of wavevectors play an important role in the physics at 
low carrier concentrations. 

In this paper, we introduce a simple one-dimensional 
quantum model which displays a multiplicity of ground 
states. Despite its simplicity, it exhibits (i) gapped states 
with commensurate density-wave order (with periods of 
2 and 3 lattice spacings), (ii) gapless regions with "float- 
ing" incommensurate, quasi-long-range density-wave cor- 
relations, and (Hi) gapped states which preserve transla- 
tional symmetry. A special case of our model appeared 
in a recent study [l| of atoms trapped in optical lattices 
2], and so an experimental study of the phase diagram 
presented here may be feasible. More generally, we of- 
fer our model as a simple laboratory, with many exactly 
solvable cases, for the interplay of density- wave orders 
with multiple periods in quantum systems. 

Our model is expressed in terms of the bose operator 
dj, which annihilates a boson on site j, and the boson 
number operator 



d)d 3 . 



(1) 



The "hard" boson condition allows no more than 1 boson 
on any pair of nearest-neighbor sites, and hence all states 
obey the constraints 



Uj < 1 



TljTlj+l = 0. 



(2) 



In the study of Mott insulators in optical lattices 0, 
the dj boson represents a dipole excitation, consisting 
of a particle-hole pair bound on nearest neighbor sites of 
the optical lattice. This microscopic dipole interpretation 



will not be crucial to our analysis here, and so we will 
refer to dj simply as a boson. 

The boson Hamiltonian we study is 



H = 



— w ( dj + dj, J + U rij + VnjUj+2 (3) 



Note that the total number of bosons is not conserved, 
and it is possible to create and annihilate bosons out of 
the vacuum. This is natural in the dipole interpretation 
of the boson, as a particle-hole pair can be created or 
annihilated from the background Mott insulator. There 
is also no explicit boson hopping term; as was shown in 
Ref. Ill boson hopping is generated by the combination 
of the constraints in and single-site terms already in 
7i, and so it is not necessary to include an explicit hop- 
ping. U is a chemical potential for the bosons, while V is 
a "nearest" -neighbor interaction, "nearest" meaning two 
sites apart, the closest two bosons can come. One can of 
course rescale out one of the couplings to obtain a two- 
parameter Hamiltonian, but it will be convenient to kee 
all three. The case V = of 7i was studied in Ref. 
we have also streamlined the earlier notation of the cou- 
pling constants to a form suitable for our analysis here. 
Without the constraints @, the Hamiltonian would be 
trivially solvable. With them, its analysis becomes quite 
intricate. 

The ground state of our Hamiltonian can exhibit 
several kinds of order, depending on the couplings. The 
Hamiltonian will favor having bosons on every other site 
if we have an attractive "nearest" -neighbor interaction, 
or a chemical potential favoring the creation of parti- 
cles. We will show that this leads to a regime of Ising- 
type order, with a spontaneously broken Z2 symmetry, 
translation by one site. If the chemical potential still 
favors creating particles, but there is a repulsive "near- 
est" neighbor potential, then the ground state will favor 
having a particle on every third site. This sort of order 
breaks a Z3 symmetry, translation by one or two sites. 



2 



When the the two kinds of ordered states have nearly the 
same energies, we will show that there exists an incom- 
mensurate phase. In the incommensurate phase, bosons 
appear on every other or every third site. 

The one-dimensional quantum Hamiltonian looks 
unusual because of the single-particle creation and an- 
nihilation operators. However, it in fact has already 
appeared in a very different context: it arises from 
taking the (Euclidean) time continuum limit of a two- 
dimensional classical statistical-mechanical model, the 
hard-square model [3J. This model describes the statis- 
tical mechanics of square tiles placed on the sites of a 
square lattice. Each square tile is rotated by 45° from 
the principal axes of the square lattice, and the area of 
each tile is twice the area of a single plaquette of the 
lattice. Tiles are not allowed to overlap, so this means 
nearest-neighbor sites cannot both be occupied: putting 
a tile on every other site of the square lattice covers all 
of space. 

It is easy to see that the Hilbert space of our one- 
dimensional quantum theory is identical to the space of 
states along a line of the hard-square model. In the quan- 
tum theory, the Hilbert space consists of bosons which 
are restricted one to a site and forbidden to be nearest 
neighbors. In the two-dimensional classical theory, the 
squares appear at most one to a site, and are forbidden 
to be nearest neighbors. It is less obvious that one can 
obtain the quantum Hamiltonian @ by taking a limit of 
the classical transfer matrix. We show in the Appendix 
how to do this, by taking the (Euclidean) time direction 
to be along the diagonals of the square lattice. Roughly 
speaking, the chemical potential U and interaction V cor- 
respond to a chemical potential and interaction strength 
for the squares. The precise relations are derived in the 
Appendix. 

The phase diagram of the hard-square model has been 
studied in Ref. yj. We will explain in section II how 
these results apply to the ground state of our quantum 
Hamiltonian. An extremely useful result is that the hard- 
square model is integrable for some values of the Boltz- 
mann weights |3( ■ In terms of the Hamiltonian @ , this 
amounts to two lines in our two-parameter space of cou- 
plings: 



uv + v 2 



(4) 



There is a single critical point along each of these lines; 
they are in the universality classes of the Ising tricritical 
point and of the three-state Potts model critical point. 
We will show how these critical points describe transi- 
tions to and from the Z2 and Z3 ordered phases. 

In section |H] we discuss the different regions of phase 
space. In section UTTl we discuss in depth the phase with 
order of period 3, the phase with spontaneously-broken 
Z3 symmetry. We present numerical results to support 
the analytic arguments. In sections llVl and IVl we discuss 
the incommensurate "floating" phase in depth. In section 
II VI we derive an effective Hamiltonian for this region, and 
prove that an incommensurate phase exists by bosonizing 
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FIG. 1: Ground-state phase diagram of Tt as a function of 
U/w and V/w. Schematic representations of the ground states 
with density waves of period 2 and 3 are shown. The dotted 
lines indicate the positions of the integrable lines in Q. The 
critical points along these lines I6I8H are labeled M2,3. The 
extent of the incommensurate region has been greatly exag- 
gerated for clarity: it is drawn to scale in Fig[S] There could 
also be an incommensurate region adjacent to the chiral tran- 
sition line about M3, but we have no direct evidence for this 
possibility. 



it. In section we calculate the size of this phase by 
using the Bcthe ansatz to compute the effective Luttinger 
parameter exactly. This enables us to find the region 
in which vortex-type perturbations are irrelevant and do 
not destroy the incommensurability. Unfortunately, this 
region occupies a fairly small region of parameter space, 
so it is not possible to see this result numerically. 



II. PHASE DIAGRAM 

Our main results are summarized in the phase diagram 
shown in Fig^ We discuss the phases and phase tran- 
sitions in this figure by considering a number of limiting 
cases. 



A. Z2 order 

First we consider the case V = 0, which was stud- 
ied earlier in Ref. QJ Although the Hamiltonian looks 
simple in this limit, the prohibition of nearest-neighbor 
bosons means that this case is not solvable except as 
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U/w — ► ±00. For large positive U/w we obtain a fea- 
tureless ground state with a small density of bosons and 
which breaks no translational symmetry, as is easily seen 
from TL. This is a disordered, or liquid, phase. Similarly, 
for U/w large and negative, the energy is minimized by 
states with a maximal density of bosons: there are 2 de- 
generate states of this type, each with every second site 
occupied; in other words, there is a density wave of period 
2, and translational symmetry is spontaneously broken. 
It was shown numerically in Ref . Q] that there is a second- 
order Ising critical point which separates these phases at 
U/w = -1.308... . 

One can also find the location of the transition from 
Ising order to disorder easily in the limit V — > — 00 and 
U rs —V. With \U\ and \V\ both large, we can ignore 
the w term, and then 7i becomes a model of classical par- 
ticles. The empty state with no particles has energy 0, 
while the maximally occupied states with density waves 
of period 2 have energy per lattice site of (U + V)/2. 
There is therefore a first-order transition between these 
states at U — — V. Actually, it is quite easy to find the 
position of this first-order transition to order (w/V) 2 . 
Using standard second-order perturbation theory in w, 
the energy per site of the empty state is —w 2 /U. Simi- 
larly, the energy per site of the maximally occupied state 
is (U + V) /2 + w 2 / (2(U + 2V)) . Equating these two ener- 
gies we estimate the position of the first-order transition 
to be 



U 
V 



-1 



w 
V 2 



(5) 



The integrability of the model along the lines (£IJ lets us 
understand how this transition occurs at finite V and w. 
In fact, the condition J5j coincides precisely with one of 
the integrable lines in ffl. Indeed, the arguments of Huse 
and the exact results [3, Q imply that integrable line is 
precisely along the co-existence line of these two phases, 
and hence the first-order boundary is given exactly by 
(JSJ), with no corrections at higher order in w/V. 

At the value 



V 

w 



V5 + 1 



5/2 



(0) 



there is a critical point along the integrable line. Thus 
the first-order transition line must terminate. Huse ar- 
gued that this integrable quantum critical point was a 
tricritical point separating a first-order boundary from a 
line of second-order Ising transitions. Indeed, the com- 
putation of exact critical exponents at this point 0, Q 
shows that this critical point in the hard-square model 
is in the same universality class as the tricritical point 
in the Ising model. This tricritical point in the contin- 
uum limit is described by the minimal conformal field 
theory with central charge c = 7/10 @. Therefore, the 
line of first-order transitions determined by JSJ) is present 
for V/w < -((y/E+ l)/2) 5 / 2 , where it terminates at the 
Ising tricritical point Mi shown in Fig^ There is a line 



of second-order Ising transitions on the other side of M2, 
but this is not on the integrable line this second-order 
line clearly passes through the point © discussed above. 
The exact computations show that the integrable line for 
> V/w > -((V5+ l)/2) 5 / 2 is in the Z2 ordered phase 
0, . This is consistent with the finite value 10 of U c 
for V = 0; at V = the integrable line is at U = —00. 



B. Z3 order 

Consider the case V — > 00, \U\ -C V. The large posi- 
tive value of V now prohibits bosons separated by 2 lat- 
tice spacings, and hence we effectively have the constraint 



njrtj+2 = 



(7) 



along with the constraints (J5J. The physics in this re- 
gion can be understood using reasoning parallel to that 
in Section III Al As before, for large positive U/w, we 
obtain a featureless ground state with a small density of 
bosons, which breaks no lattice symmetry. For U/w large 
and negative, the ground state has a maximal density of 
bosons: with the constraints J3J) and (JJJ there are now 3 
degenerate states of this type, each with every third site 
occupied. Here we thus have a density wave of period 3. 

To study this phase with spontaneously-broken Z3 
symmetry at finite V, we can use results for the other 
integrable line, where V/w > in JSJ). Along this line 
there is a quantum critical point at |3J 



V 

w 



V5 + 1 



5/2 



(8) 



This critical point is in the universality class of the three- 
state Potts model 3, 5]. In the continuum limit this is 
described by the minimal conformal field theory with cen- 
tral charge 4/5 @. For w/V smaller than this value and 
along this line, the exact results indeed show that the 
Z3 translation symmetry is broken. For w/V larger than 
this value and along this line, there is no order. Thus 
along the integrable line, the transition is in the univer- 
sality class of the usual order-disorder transition in the 
three-state Potts model; similar Z3 quantum criticality 
occurs in an XXZ chain in a staggered field Q- However, 
as opposed to the Ising case, the transition away from 
this integrable line is not of the usual type £| : it is re- 
lated that of the chiral clock model @. We will discuss 
the nature of the transition between these states as a 
function of U /w in Section IIIII including the possibility 
of additional intermediate states with incommensurate 
periods. 



C. Competing orders 

The most complicated and interesting limit is V — > 00, 
U ~ — 3V. where there is competition between density- 
wave orders of period 2 and period 3. Initially, as in 
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Section fll Al we anticipate that with \V/w\ and \U/w\ 
both large, we can neglect quantum effects and consider 
a model of classical particles. The period-2 and period-3 
density- wave states then have energies per site (U + V)/2 
and U/3 respectively, so that the ground states become 
degenerate at U = —3V. In fact, at U = —3V in this 
limit, there are many degenerate states: any state where 
successive particles are separated by either 2 or 3 lattice 
spacings is within the ground state manifold. This results 
in an extensive classical entropy at U — — iV . This ex- 
tensive degeneracy is lifted only by quantum effects, and 
perturbation theory in w is therefore highly non-trivial. 
The analysis of quantum effects in this regime will be 
presented in Sections IIVI and E] We will establish that 
there is an intermediate gapless phase with incommen- 
surate spin correlations. We will show as well that the 
gapped state with no broken symmetry also intervenes 
between the period-2 density-wave state and the incom- 
mensurate phase, as shown in Fig. ^ 

We will argue in Section lUfl that there also should be 
an incommensurate phase near the Potts critical point 
M 3 . The simplest and most plausible possibility is that 
this phase and the one for — U « 3V 1 are the same. 
We have drawn the phase diagram in this fashion. 



III. DENSITY- WAVE ORDER OF PERIOD 3 

The density-wave state of period 3 appears when V 
is large and positive, so that bosons with separations of 
both 1 and 2 lattice spacings are suppressed. This section 
and the next will discuss the manner in which quantum 
fluctuations destroy this period-3 order. 

It is useful to begin with the limit V — > oo, noted 
above in Section lllfjl At first glance, there appear to 
be only two distinct possibilities for the ground states: a 
featureless, low density state at large positive U/w, and 
a density- wave state of period 3 at large negative —U/w. 
We might also anticipate that these states are separated 
by a second-order critical point in the universality class of 
the 3-state Potts model. However, the situation is more 
subtle Hi- 
lt is useful to introduce order parameters, H! p , which 
characterize density-wave orders of period p (in general, 
p need not be an integer): 

3 

For p = 2, the Ising order parameter ^ is real. However, 
for the p = 3 case we consider in this section, the Potts 
order parameter $3 is complex. Clearly, the state with 
density wave order of period 3 has (^3) ^ 0, while the 
featureless, low density state has (^3) = 0. 

Using the usual symmetry criteria, we can write down 
a continuum quantum field theory describe the onset of 
Potts order. The action should be invariant under trans- 
lations, under which ^ 3 — > e l2n ^ 3 ^3, with I integer. A 



crucial point, noted by Huse and Fisher |9j in a different 
context, is the behavior under spatial inversion: 

X^-X , * 3 ^*3 (10) 

where x is the continuum spatial co-ordinate. Finally, we 
note that 'J 3 is invariant under time reversal. These sym- 
metry constraints lead to the following proposed effective 
action for the quantum field theory (r is imaginary time) : 

S 3 = J dxdT[iay* 3 d x y 3 + c.c + \d T y 3 \ 2 +v 2 \d x V 3 \ 2 

+ r|* 3 | 2 +i>^ + c.c. + ...] (II) 

The critical point associated with this field theory occurs 
at a = r = 0, and is in the same universality class as 
the critical point in the three-state Potts model. This is 
described by a conformal field theory of central charge 
4/5, and the dimensions of all the operators are known. 

The exact results show that the integrable critical 
point M 3 indeed is in the universality class of the three- 
state Potts model 0,0 ■ It separates a state with density- 
wave order of period 3 from the gapped translationally- 
invariant state 0, |(| . Just because the critical point is 
in the same universality class as the three-state Potts 
model does not mean that all the physics in the region of 
the critical point is the same. The key difference is the 
presence here of the linear spatial derivative proportional 
to a, which is clearly allowed under the symmetries IIOII . 
This term is chiral, in that it breaks rotational symmetry. 
The exact results show no evidence of chiral behavior, so 
the effective theory describing the integrable line presum- 
ably corresponds to a = 0. Varying r at a — therefore 
moves one along the integrable line, and describes the 
usual physics of the three-state Potts model: a second- 
order transition between Z3 order and disorder. 

When one moves off the integrable line (taking a ^ in 
the effective theory), the theory in the universality class 
of the chiral clock model 0. The chiral perturbation, 
a, is known to be relevant at the ordinary Potts critical 
point with scaling dimension 1/5; the perturbing op- 
erator has dimensions 9/5 and breaks Lorentz invariance. 
Thus for a^0, the quantum phase transition associated 
with the vanishing of density-wave order of period 3 is 
in general not in the same universality class as the usual 
order/disorder transition in the ordinary 3-state Potts 
model. The integrable point M 3 is an isolated multicrit- 
ical point; all other points on the phase boundary of the 
period-3 state must have a^0, and display the physics 
of the chiral clock model. Note the contrast with the cor- 
responding transition associated with density-wave order 
of period 2: this transition is in the Ising universality 
class because 'J 2 is real, and then the term proportional 
to a is a vanishing total derivative. 

The chiral clock model is generally believed to exhibit 
an incommensurate phase 0, Il0j| . In the region near 
the multicritical point M 3 , this occurs when the chiral 
perturbation (a 7^ 0) dominates over the thermal per- 
turbation (r /: 0). This region is therefore quite narrow, 
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because a, of dimension 1/5, is much less relevant than 
r, which is of dimension 6/5. Thus r must effectively 
vanish for the effects of the chiral perturbation to be felt. 
In the next section, we will show that the incommensu- 
rate phase does have a finite (albeit small) width far from 
the multicritical point, where \U\ and V are large, with 
U =~ V/3. There is no indication that any more new 
physics intervenes in between this region and the incom- 
mensurate region near M3 , so we presume that these two 
phases are the same, as we have indicated in Fig. ^ 

In the remainder of this section we will describe our nu- 
merical study of the physics in the limit V — > 00, which 
corresponds to being on the portion of the phase bound- 
ary of the period-3 state above the multicritical point M3 
in Fig^ Our numerical results here are inconclusive, as 
we see little deviation from the pure 3-state Potts physics 
for the system sizes examined. We will study the portion 
of the phase boundary below M3 in Section llVl and there 
we shall definitively establish the existence of a gapless 
incommensurate phase, driven by the presence of the chi- 
ral perturbation; the ordering wavevector, K = 2ir/p, of 
this incommensurate phase obeys K > 2ir/3. A simple 
fluctuation analysis of £3 indicates that K — 2n/3 is pro- 
portional to a. As a vanishes at M3 where K = 27r/3, 
these results suggest the conjecture that there is a small 
intermediate, gapless, incommensurate phase also above 
the point M3, but with a change in sign of a leading to 
an ordering wavevector K < 27r/3. Our numerical stud- 
ies below, however, do not show any specific evidence 
in support of this conjecture. This is possibly because 
the scaling dimension of the chiral perturbation is quite 
small, and very large system sizes are needed before its 
effects are perceptible. Furthermore, we shall find in the 
study of the portion of the phase boundary below M3 in 
Section IIVI that the region of incommensurate spin cor- 
relations is extremely small, and it is likely that a similar 
feature holds in the above M3 region being studied here. 

The numerical analysis for large V is simplified by sim- 
ply eliminating all states with an energy determined by 
V. This is equivalent to extending the constraints (01 to 



Uj < 1 



IjTlj+l 







0. 



and working with the simplified Hamiltonian 



H s = 



[- w (dj+dfj+Unj 



(12) 



(13) 



The exact diagonalization of 7i p can then be carried out 
for system sizes N < 21. 

First, we searched for the preferred values of p in a 
system of size N = 21. We applied a small external 
cosine potential on the bosons of magnitude 10~ 4 and 
wavevector 2w/p, and measured the resulting values of 
(iffp) as a function of U/w. The results are shown in 
Figs |2] and Fig. shows that value of $3 is much 
larger than the response at another nearby wavevector. 
Figure |3 shows similar data, but with a greatly expanded 
vertical scale. The response at periods p ^ 3 is always 




FIG. 2: Density-wave order, (9 P ), induced by an external 
potential at wavevector 2tt/p and magnitude 10 -4 . There 
appears to be onset of spontaneous density-wave order with 
period p at the smaller values of U/w. 
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FIG. 3: As in Fig. [5] but with a greatly expanded vertical 
scale. All of the periods shown above show no sign of spon- 
taneous order. 



much smaller, and so we see no sign of spontaneous or- 
der at any other wavevectors. Interestingly, the largest 
secondary response is at the wavevector 167r/21. This 
is contrary to our expectations earlier based on the in- 
fluence of the chiral perturbation in (|ll|l . where we an- 
ticipated incommensurate order at wavevectors smaller 
than 27r/3. However, given the very small susceptibility, 
and the small system sizes considered we are not able to 
interpret any of this data in support of an incommen- 
surate phase. Any such phase, if present, must be over 
an extremely small parameter regime, as is found in Sec- 
tionllVl 
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FIG. 4: Plot for the energy gap A as a function of U for 
different system sizes. The data for different system sizes 
cross at U c /w — —1.852. 
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FIG. 5: Scaling plot for the energy gap A. The data for dif- 
ferent system sizes collapse for v — 0.83 and U c /w — —1.852. 



Ignoring the presence of a possible intermediate in- 
commensurate phase, we now see if our data can be in- 
terpreted in the context of a direct transition from the 
gapped density-wave state of period 3, to a gapped "liq- 
uid" state which does not break any translational sym- 
metries. We expect such a transition to be in the Potts 
universality class. Indeed, we now show that our data 
are in complete accord with such an assumption, with 
excellent agreement to the expected values of the critical 
exponents. This is again consistent with an as yet unde- 
tected influence of the chiral perturbation in the present 
analysis (the chiral perturbation is more clearly detected 
in Section HVJl. 

First, let us determine the position of a presumed crit- 
ical point from the period-3 (Potts) density-wave state to 
a disordered state using finite size scaling analysis. Near 
the quantum critical point, the energy gap A is expected 
to obey the scaling relation 



A = N~ z <!> 



(N^"(U-uSj (14) 



where $ is an universal scaling function, z is the dynamic 
critical exponent, and v is the correlation length expo- 
nent . Fig. 0]shows a plot of N A / w vs UJ w which exhibits 
a clear crossing point at 



0.36 



0.32 



0.28 



0.24 



OT 0.20 



0.16 



0.12 




-1.92 -1.90 -1.88 -1.86 -1.84 -1.82 -1.80 -1.78 
U/W 

FIG. 6: Plot of the equal-time structure factor S2n/3 for dif- 
ferent system sizes. The data for different system sizes cross 
at U c /w = -1.852 for r\ = 0.27 



Potts order parameter '5 3 (Eq. EJl given by 



U c 

— = -1.852. 

w 



(15) 



From Eq. El we thus conclude that z = 1. Fig shows 
a plot of NA/w vs iV 12 ([7 - U c )/w. We find that the 
data for all N collapses for U c /w = —1.852, giving v — 
0.833. The numerical values of v and z obtained here are 
in excellent agreement with the known analytical values 
z = 1 and v = 5/6 for the three-state Potts model [l7| . 

As a final test of the three-state Potts universality 
class, we calculate the equal-time structure factor of the 



2tt/3 



N 



(16) 



The structure factor is expected to scale as N 2 ~ z ~ v near 
the critical point. A plot of «S 2w /3/A'" C1 - 73 vs U/w, as 
shown in Fig. |H| shows excellent crossing at U c /w = 
— 1.852. From this, knowing z = 1, we identify 77 = 0.27, 
which is again i n g ood agreement with the known exact 
value 77 = 4/15 |l7j . 
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IV. COMPETITION BETWEEN 
DENSITY- WAVE ORDERS OF PERIODS 2 AND 
3: THE INCOMMENSURATE STATE 

We noted in Section III CI that in the classical limit 
(w = 0), the period-2 and period-3 states become de- 
generate for U = —3V. Here we shall study the quan- 
tum physics near this degeneracy point. As will become 
clear, the physics is quite complicated, but we are able 
to obtain a complete picture by a combination of pertur- 
bation theory, bosonization, and Bethe ansatz methods. 
Our primary result is that the gapped period-2 and pe- 
riod 3 are separated by 2 intermediate phases, and all 
these phases are separated by 3 second-order quantum 
phase transitions (see Fig.^J. Above the period-2 phase 
is the gapped, translationally invariant state, accessed 
by a conventional Ising transition. Above this is a gap- 
less, incommensurate phase reached across a Kosterlitz- 
Thouless transition. Finally, we reach the gapped period- 
3 phase via a Pokrovsky-Talapov transition 11]. 

We will establish these results by a careful study of the 
region of the phase diagram where 



O 



o o 



o o 



o 



o 



0<w,\U + 3V\ < \U\,V 



(17) 



It will become clear from our analysis that the non-trivial 
region of the phase diagram is where 



\U + 3V\~—. 



(18) 



So we introduce the new dimensionless coupling a, and 
parameterize 



U 



-3V 



1 + 


/ w\ 2 









(19) 



The above conditions on the interesting regime can now 
be restated as follows 



< w < V 



1 



(20) 



We shall follow the physics as a is tuned from large neg- 
ative to large positive values, while maintaining j!7jl. 



A. The effective Hamiltonian 

As we mentioned in Section III CI under the condition 
l|2Up. the important states are those in which successive 
dj bosons are separated by either 2 or 3 lattice spacings. 
We denote such a state by the sequence of lattice spacings 
e.g. | . . . 23323233 . . .); see Fig.0 We wish to develop an 
operator formalism for describing such states. We shall 
introduce two distinct formalisms below, but it should al- 
ways be kept in mind that both describe identical Hilbcrt 
spaces with identical spectra. 

In the first approach, we describe the states in terms 
of 'defects' or 'domain walls', or 'kinks' in the ordered 
period-2 density-wave states. This description focuses 



FIG. 7: The state | . . . 23322 . . .). Filled circles are dj bosons, 
and empty circles have no bosons. 



on the 3's in a background of 2's. We view each 3 as 
a new boson, created by the operator t\, located at the 
center of the 3 defect. Note that the I sites lie at the 
midpoints of the links of the j sites of the original model 
H. Thus the state | . . . 2223222 . . .) represents a single 
t boson. An arbitrary state in the low-energy subspace 
can be described by a configuration of t bosons. The 
resulting t boson states obey some constraints: there can 
only be at most one t boson per site, and successive t 
bosons must be separated by 3, 5, 7, 9, 11... sites. 

Second-order perturbation theory also allows us to ob- 
tain the effective Hamiltonian for the t bosons to order 
(w/V) 2 . It is useful to first tabulate the energies of dif- 
ferent configurations of the original d bosons in different 
environments, to order {w/V) 2 . For a d boson between 
22 bonds we have 

2 

E 22 = U + V + 



U + 2V 

Similarly for a d boson between 33 bonds we have 

■" 2 

E33 



U + TT' 



Finally for a d boson between 23 bonds we have 

,r 2 

E23 



ur 



V 

U+ 2 ' U+V 



From these expressions, we can obtain the energy of a t 
boson well-separated from all other t bosons: 



E t — 2E23 — -E22 



V 



3(1 + a) (w 



Similarly, we obtain the interaction energy of two adja- 
cent t bosons 



IE' 



•13 



= V 



+ E33 

1 iw 
3 \V 



4E 9 



O 



2Ek 
3" 



(v) 



Finally, moving on to the t boson hopping amplitude, a 
similar perturbation theory shows that it is w 2 /{U + V). 
Collecting all these results, we can now write down the 
effective Hamiltonian for the t bosons, valid in the regime 
(|20|l . On a system of iV sites the Hamiltonian Ti. in (J3J) 
is effectively 

Ht = 
V 2 

I 

2 



(1 



+ (|0 E -^(4 +2 ^+4^+0 



(21) 



3, . + 1 j j 
-(1 +a)t £ tt — -t^ 3 U + st e t£ 
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We reiterate that Ht must be solved under the constraints 
that there can only be at most one t boson per site, and 
successive t bosons must be separated by 3, 5, 7, 9, 11 . . . 
sites. To this order in w/V the Hamiltonian conserves 
the number oft bosons. Except for the latter constraint, 
this Hamiltonian is very much like that of an XXZ mag- 
net. The hopping term in (1221 is akin to spin exchange, 
the chemical potential is like an effective magnetic field 
h e ff = 3(l + er)/2, and the interaction is akin to the S Z S Z 
interaction between neighboring sites. The analogous 
(ferromagnetic) XXZ model has anisotropy J z /J x = 1/3, 
so it is in a gapless regime. We shall show in sect. [V] 
that this model is as well, for appropriate values of the 
effective magnetic field. 

Before turning to analysis of Ht, it is useful to intro- 
duce a complementary analysis of the physics valid under 
the conditions (|20|l . Now, we consider defects or kinks 
between ordered states with period 3. So we treat the 
2's as bosons moving in background 3's. Such bosons are 
created by the operator pj , and note that this boson re- 
sides on the sites of the original lattice. The energy of 
an isolated p boson is 



E n — 2E- 



23 



3^33 



V 



Also, the interaction energy of two adjacent p bosons is 



E 



p,int 



= AE 



23 



Eon — AE. 



33 



3E n 



V 



•-fSV 

3 \VJ 



Finally, the hopping matrix element of the p bosons is 
the same as that of the t bosons. Collecting these re- 
sults yields the following effective description of H in the 
regime (|2"U)l : 



Hp 



(1 + 9cr) (W^~ 

,V 



1 



2 1^+3^" 



-PjPj+3 



(22) 



The constraints on the p bosons are that each site can 
have at most one p boson, and successive p bosons can 
only be separated by intervals of 2, 5, 8, 11, 14, . . .. Note 
that except for the constraints, this description also re- 
sembles a (ferromagnetic) XXZ model with the same 
anisotropy A = J z /J x = 1/3. 

It is important to note that the Hilbert space and spec- 
tra of Ht in l(2^|) and H p in lt2"5|) should be identical, 
although this is by no means obvious. Some simple con- 
sistency checks can however easily be performed: the vac- 
uum state of Ht with no t bosons should have the same 
energy as the state of H p with the maximal number of 



z 3 

ordered 

oo«oo«oo 



incommensurate 



disordered 



-1 



z 2 

ordered 



a 



n ► 





FIG. 8: Ground states of Ti, Tit, or TL P as a function of a, 
under the conditions Il2()> . The dimensionless coupling a is de- 
fined in l|19fl : so the horizontal axis above corresponds to mov- 
ing downward along the extreme left of Fig0 The transition 
from the Z2 ordered state to the disordered state at a = —1/3 
is in the Ising universality class, and has dynamic exponent 
z — 1. The Z3 ordered state first undergoes az = 2 transi- 
tion to an incommensurate phase in the Pokrovsky-Talapov 
universality class at a — —13/9, which is then followed by a 
z — 1 Kosterlitz-Thouless transition at a = —1.422... (this 
last number is determined from the Bethe ansatz analysis in 
Section EJ- 



p bosons. The reader can easily check that this require- 
ment, and its converse, do indeed hold. More generally, 
we can easily see from the manner in which we defined 
the states that if a state has N t t bosons under Ht , then 
the same state has N p p bosons under H p where 



2N P + 3N t = N, 



(23) 



and N is assumed to be a multiple of 6. 

Some simple considerations now allow us to deduce 
some important properties of Ht and H p . Consider first 
the case <r> 1. A glance at Ht shows that each t bo- 
son costs a large positive energy, and so the ground state 
is the t boson vacuum (which also happens to be the 
state with the maximal number of p bosons). This vac- 
uum state is clearly the state with density state of period 
2. The lowest excited state consists of a t boson excita- 
tion (a kink), at zero momentum, which has an energy 
(w 2 /V)(l + 3ct)/2 above the ground state. So we con- 
clude that the ground state remains a density wave of 
period 2 provided 



a > 



period-2 density-wave order 



(24) 



This conclusion is illustrated in Fig. [SI which displays 
the final phase diagram emerging from the analysis of 
this section. 

A similar argument can be applied to determine the 
stability of the period-3 density-wave state. From H p , 
this is stable for cr <S — 1. The energy of a single p boson 
excited state at zero momentum is — (w 2 /V)(13 + 9<t)/2, 
and so the ground state remains a density wave of period 
2 provided 



cr < -- 



13 



period-3 density-wave order (25) 



This is also shown in Fig. [S] 
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B. Stability of the incommensurate phase 

It now remains to understand the physics for —13/9 < 
a < —1/3. The above considerations suggest that this 
intermediate region has a finite density of both t and p 
bosons obeying J33J). If this is a compressible phase, then 
we expect the ground state to be a Luttinger liquid, and 
the corresponding correlations of the density of the un- 
derlying d bosons to be incommensurate. Furthermore, 
the transition a = —13/9 involving the onset of a non- 
zero density of the p bosons would be a transition of the 
Pokrovsky-Talapov type 0] > with dynamic critical expo- 
nent z = 2 12], and an effective free fermion description 
near the critical point. Similar considerations can be ex- 
pected to apply to the transition at a — —1/3, associated 
with the onset of a finite density of t bosons. 

However, the validity of such a conclusion requires the 
absence of additional perturbations, at higher order in 
w/V ', which could disrupt the stability of the Luttinger 
liquid phase. Of particular importance here is the pos- 
sible generation of 'vortex operators' [l3l IbH flU ] which 
violate the conservation of the total number of t or p 
bosons obeyed by the terms so far included in TLt, P - A 
simple analysis of the processes at higher orders in w/V 
in terms of the b bosons shows that such processes do 
indeed appear: the domain walls in the period-2 density- 
wave state, the t bosons, can be annihilated or created 
in pairs. Similarly for the period-3 state, there are pro- 
cesses in which 3 p bosons can be annihilated or created. 
We will examine the consequences of such processes by a 
combination of bosonization and Bethe ansatz methods. 
We will find that the vortex perturbations are indeed ir- 
relevant at the a = —13/9 critical point, so that there 
is a direct z — 2 transition from the gapped period-3 
density-wave state, to a stable, gapless incommensurate 
phase. However, the vortex perturbations are found to 
be relevant at the a = —1/3 critical point, implying that 
the gapless incommensurate phase does not extend all the 
way up to the period-2 density- wave state (see Fig- - 

Here we study the stability of the Luttinger liquid 
phase of Tlt, P postulated above within the framework of a 
long-wavelength bosonization analysis. This will allow us 
to place general conditions on the Luttinger parameters 
required for stability. This analysis will allow us to de- 
termine the values of the Luttinger parameters in some 
important limiting cases, and show there indeed exists 
an incommensurate phase stable to perturbation. The 
analysis does not allow us to determine the Luttinger 
parameter in general — this we will do in the following 
section by a Bethe ansatz analysis. 

First, let us carry out the analysis using H. t , and let us 
assume the value of a is such that we are in a compressible 
Luttinger liquid region. We follow the notation of Ref.lili 
and introduce the continuum field 9 t with the action 

Here v is a velocity and K t is a Luttinger parameter, 



normalized so that free fermions have free fermions have 
K = 1. In terms of this field 

t~e- i6 K (27) 

Now the vortex operator V ~ t 2 ~ e ~ 2lB t j s seen ^ have 
scaling dimension 1/Kt under the action St, and hence 
the vortices are irrelevant, and the compressible phase is 
stable, as long as 

Kt<\- (28) 

We will determine the value of K t as a function of a 
in the following subsection, but here we note a simple 
where K t can be determined exactly. This is the limit 
where the density of t bosons becomes vanishingly small 
as we approach the period-2 density-wave state with a /* 
— 1/3. In this low density limit near a z = 2 quantum 
transitions, it can be shown that the t particles behave 
like free fermions 0] 1 an d so we must have K t = 1. The 
same argument can also be used to determine v and so 
we have the important result 

v = {w 2 /V)V~l - 3cr ; K t = l as a/ -1/3. (29) 

Notice that Kt does not obey l|28(l , and hence vortices are 
always relevant at the boundary of the period-2 state, as 
we claimed above. As these vortices correspond to cre- 
ation of pairs of the free fermions, the resulting fermion 
model is easily seen to be equivalent to the fermionized 
Ising model. This implies that we have an Ising tran- 
sition at a = —1/3 directly from the gapped period 2 
ordered state to a disordered gapped state that does not 
break any symmetries, as shown in Fig|SJ 

Next, a complementary bosonization analysis can be 
applied to H p . We express the low energy physics using 
a Luttinger field 9 p with the action 

Sp = / dxdT [ (9r0 f )2 + v2 ( d * p) 2 ] ( 3 °) 

As TLt and TL p have the same spectrum, the velocity v 
should be the same as that in (|26[l . However, the Lut- 
tinger parameter K p is now different because the inter- 
pretation of the field 8 p differs from (|2*7Jl : 

p ~ e- ie " (31) 

Now the vortex operator ~ p 3 has scaling dimension 
9/(4fT p ), and hence the vortices are irrelevant, and the 
compressible phase is stable, as long as 

K P < I (32) 

Again, determination of the values of K p and v requires a 
Bethe ansatz analysis, but exact results can be obtained 
in the limit a \ —13/9, where the density of p bosons be- 
comes vanishingly small. The same arguments as applied 
above for the t bosons now show that 

v = (w 2 /V) v / 2{13 + 9a)/9 ; K p = 1 as a \ -13/9. 

(33) 
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Now K p does obey the stability condition (|32(l . and hence 
vortices are always irrelevant at the boundary of the tran- 
sition to the period-3 density-wave state. So, as claimed 
earlier, we have established the existence of a gapless, 
incommensurate state in this region. At some interme- 
diate value of (j, the vortices will become relevant and 
will drive a Kosterlitz-Thouless transition to a gapped, 
disordered phase. 

We conclude this subsection by noting that the rela- 
tionship H23[l implies a simple, general relationship be- 
tween the Luttinger actions in (|26[l and (|3U[) . Indeed, 
the standard connection between particle density and the 
9 field in bosonization theory [0|, combined with (|23fl . 
leads to the condition 



(34) 



This result is also consistent with our expressions for the 
vortex operator V noted below l|27|l and l|31|l . Using 134(1 
in 126(1 and l)3()[l. we also obtain the exact relationship 



(35) 



which shows, not surprisingly, that the conditions 1(28(1 
and 1(32(1 are the same. 



V. THE SIZE OF THE INCOMMENSURATE 
PHASE 

To establish the boundaries of the incommensurate 
phase, we will show that the effective Hamiltonian 1(22(1 
(or equivalently, can be solved using the Bethe 

ansatz. We will then use this to determine the Luttinger 
parameter anywhere in the incommensurate phase. The 
computation is a generalization of that done in |20| to 
cases with a "hard-core" constraint. Bethe ansatz anal- 
yses have been done in several related cases, the XXZ 
model with a hard core (not allowing down spins to be 
adjacent') 1 18|. or a model of hard-core fermions arising in 
a supersymmetric chain |19| . The details of the computa- 
tion of the Luttinger parameter have not been presented, 
however, so we describe the calculation in generality here. 



A. The Bethe ansatz 

The t boson in the effective Hamiltonian (|22|) is a kink 
which separates the two density-wave states. We denote 
the Ising-ordered state with no kinks as |0), so that an 
eigenstate of the Hamiltonian (|37|) with N t < N/3 kinks 
is 



The ij are ordered so that 1 < i\ < %i < 



over we require that 



3,5,7, 



4 10) (36) 



. ., and more- 
the Hamilto- 



nian preserves this restriction. To simplify the analysis 



in this section, we shift away the constant term in Tit, 
multiply by a constant, and define the effective chemical 
potential h t — — 3(1 + cr)/2. This yields 



H 



E 



t i+2 ti + tgte+2j + httfti + -t e+3 ti + 3t e tg 

(37) 

This Hamiltonian has the same eigenvectors as the orig- 
inal. For simplicity, we also assume periodic boundary 
conditions. 

Bethe's ansatz is that the state <j>( N *) is of the form 



(38) 



for some numbers /ij and Ap, where j = 1, . . . , N t and 
P = (PI, P2, . . . , PN t ) is a permutation of the integers 
(1, 2, . . . , Nt). With periodic boundary conditions, we 
can construct eigenstates of the translation operator T, 
which sends Tt\t\ . . . = 4i+i4 2 +i ■ ■ •■ A state obeying 
Bethe's ansatz is an eigenstate of T if the amplitudes Ap 
are cyclically related as 



ApN t ,Pl,P2,...P(N t -l) — Mpjv t ^fi,f 2, 

If we define the bare momenta kj via /Xj 
eigenvalue e zktot of T is then given by 



■ PN t - 



(39) 



e lkj , the 



fctot = -i ha 



N t 



(40) 



The relation l|39|) then can be thought of as the quanti- 
zation of the momentum of a particle in a box. 

The fii and the amplitudes Ap are found by demanding 
that cj)^ be an eigenstate. It is simplest to illustrate this 
in the case of two kinks. The Bethe ansatz is that 

¥>{h,i2) = E ^ p MpiMpa 
p 

= A 12 tftf +A 21 ^tf 
Requiring (f>^ be an eigenstate of T means that 
A u = Aai/xf , 

with e~ lktat = /ii/i2- Operating with the Hamiltonian 
(|37|l on 0( 2 ) yields 



H<t> 



(2) 



E f 6^ 



where 
E t 

and 



-htN t - - [Mi+Ml + (Mi) 2 + (M 2 ) 



Xi = --<p(i, i + 3) + <p(i, i + l) + <p(i + 2,i + 3). 
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The Bethe equations are derived by requiring that all the 
Xi vanish, so that tp^ is an eigenstate. With the Bethe 
ansatz, this means that 



-Xj — — ^jAp fJ.pxfJ.p2 (2A^p 2 - fi P2 - fJ.pxfJ.p2) 
p 

where A = 1/3. Each of these vanishes if 
A 2 i fJ2 fj\fj\ + 1 - 2A/x| 







Ax 



fix nlnl + 1 - 2Afi\ ' 



Combining this with the periodic boundary conditions 
P?l> yields 



A21 



Ml 



Let i/j — fj?y Then using the fact that e tktot — fix^i 
gives the Bethe equation 



v (N-2)/2 e -ik ioi 



VrnVj + 1 - Vrr, 
V m Vj + 1 - V.j 



which holds for (m,j) = (1,2) and (2,1). One solves 
these two equations for vx and V2 subject to the con- 
= vxV2- Then the corresponding eigen- 



stramt e 



-2ifc t . 



state is found (up to an overall constant) by substituting 
the values of i/j into the equation for A21/A12. 

This computation can be generalized to any number 
of kinks N t . The constraint that cf>( Nt > be an eigenstate 
is basically the same as the two-kink case: one requires 
that 

= 2A(p(...,i,i + 2,...)-(p(...,i,i + l,...) 
-tp(...,i + l,i + 2,...) 

for any choice of ix, 12, ■ ■ ■ ^N t and i. The trick to make 
this vanish is to consider the permutation P', which dif- 
fers from P only in that Pm and P(m + 1) are reversed, 
i.e. P' = P1P2 . . . P(m-l)P(m+l)PmP(m+2) . ..PN t . 
Therefore is an eigenstate if for all P and m 



with 



A, 



g(a,b) 



g{vp{ m +l), v Pm) 



y/a(ab + 1 - 2Aa) 
Vb(ab + 1 - 2A6) ' 



(41) 



One can think of g as the bare S-matrix describing the 
phase shift when two kinks are interchanged. The fact 
that g{a,a) = —1 means that the wavefunction vanishes 
if two of the bare momenta kj are identical. To find the 
Vj, we impose the boundary condition (|39|l . Note that 



ApN t ,Pl...P{N t -\) Api,PiV t ,P2...P(jV t -l) 



1P1,P2,...PJV 4 



Apx, 



P2...PN t 



g{vpN t ,vpi) 



N t 



Y[g(vpN t ,vpj) 
J'=l 



The condition (|39(l must hold for all PNt, and hence all 
j. Putting this all together yields the Bethe equations 
for the fjj or vf 



N/2 _ 



Nt 



This is thus a coupled set of polynomial equations for the 
fii. Using the explicit form of g, and the expression l|40|l 
for the translation eigenvalue e lp , the Bethe equations 
simplify to 



Nt 



(N-N t )/2 ikto 



n 



VjV m + 1 — 2AzAj 



VjVm + 



1 - 2A^, 



These are a set of Nt coupled polynomial equations for 



the v.. 



1 . . . Nt- Solving these for a set of 



kj , one then finds the corresponding eigenstate (up to an 
overall normalization) by using (|41|l . This eigenstate has 
energy 



Ni 



E t = -h t N t - \ cos(2%) 



(43) 



These are precisely the Bethe equations and energy one 
obtains for the XXZ model at A = J z /J x = 1/3 with 
(N — N t )/2 sites, if one imposes twisted boundary con- 



ditions resulting in the factor of e 



— ik t . 



O- The 



twist and the different number of sites result from the 
excluded-volume effect (the restriction that t bosons 
must be 3, 5, 7, . . . sites apart), and the fact that t bosons 
hop two sites instead of one. 

In the next subsection we will use these Bethe equa- 
tions to find the size of the incommensurate phase. Since 
we saw by using bosonization that the incommensurate 
region is closer to the Z3-ordered phase than to the Z2- 
ordered phase, the density of p bosons in the incommen- 
surate phase will be much smaller than that of the t 
bosons. It is thus useful to derive the Bethe equations in 
the p boson basis. The analogous wavefunction tp must 
obey 

= 2A£(..., m + 2, ...)-£(..., M-l,...) 
-£(...,« + 3, i + 2,...) 

with, crucially, the same A = 1/3 as for the t bosons. 
Going through a virtually-identical calculation yields 



e 3ifctot _ - . , 



N„ 



1 \ 

E p = -h p N p - - ^ ["j + ^7 



(44) 



N v 



(N+N p )/3 iktat 



e ifctot = (_i)JVp-i yi 



LOiLOrn + 1 _ 2AiUj 



■1 bJjU>m + 1 - 2Aw„ 



(45) 
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where the effective chemical potential here is h p = a + 
4/9. Note that the right-hand-side of the Bethe equa- 
tions here are identical to those for the t bosons (and for 
the corresponding XXZ model); only the left-hand-side 
changes. 



useful because it puts the Bethe equations (|45|l in differ- 
ence form: taking the log of both sides gives 



{2M j + 1)tt = RNkiOj) - Afc tot + ^ L{9 3 



(47) 



B. Calculation of the Luttinger parameter 

Since a great deal of discussion of the size of the in- 
commensurate phase has appeared in the literature (see 
e.g. [l3l PiU Il5^. we feel it is worth discussing in detail 
a case where the answer can be determined exactly via 
the Bethe ansatz. In this subsection we therefore derive 
explicitly the size of the incommensurate phase by find- 
ing the range of a where the vortex operators are irrele- 
vant. We do this by computing the Luttinger parameters 
K p = 9K t /4 of the effective Hamiltonian. The dimension 
of operators can be computed directly by using either 
the asymptotics of the correlator [2l| or finite-size effects 
[22j| . but the result is identical. In this subsection, we 
will generalize Haldane's computation of the Luttinger 
parameter |2(il | to the case at hand. The computation 
we do here is more complicated because the exponent on 
the left-hand-side of (|42|l or (|45() depends on the number 
of kinks in the system. This is a consequence of the re- 
strictions on the locations of the particles: more particles 
effectively reduces the size of the system. 

The XXZ spin chain at A = 1/3 is in a gapless phase. 
The vortex operators are forbidden from occuring be- 
cause of the t/(l) symmetry of the spin chain; if one 
breaks this symmetry by allowing e.g. J x ^ J y , then 
the resulting XYZ chain is indeed gapped. One can show 
using the above Bethe ansatz computation that Hamilto- 
nians H'2'2\) and are also in a gapless phase: the slight 
difference in Bethe equations does not change this result. 
This difference, however, does complicate the computa- 
tion of the dimension of the vortex operators. It also 
means that as opposed to the XYZ model, there is a 
gapless phase with broken U(l) symmetry and irrelevant 
vortex operators. In a gapless phase, the Luttinger pa- 
rameter Kp is defined by the general relation ^(j 



where Mj is an integer, and 



K„ 



d 2 E p 



(46) 



where E p is given in l|44|) . and v is the velocity of the 



2N n 



N means 



excitations. Using the fact that 3N t 
that K p = 9K t /A, as noted above. 

To do this computation, it is convenient to rewrite 
the Bethe equations in terms of "rapidity" variables 9j. 
These are defined by the relation 

-3ik(6j) sinh(% + ij/2) 
smh{9j - vy/2) 

where 7 = cos -1 (—A); for our models, A = 1/3. This 
relation also defines A: as a function of 9, i.e. kj = k(9j). 
This change of variables from momentum to rapidity is 



L(9) = i In 



sinh(0 + vy) 
sinh((9 — ij) 



We have defined the parameters R and A to make the 
analysis apply to general models with hard-core con- 
straints. For the p bosons we have R = 1 + N p /N and 
A = 1. For the t bosons, we have R = 1 — N t /N and 
A = -1. The usual XXZ model has R = 1 and A = 0. 
The cases discussed in [Til Il9j can also be written in this 
form. 

Each kink is characterized by an different integer My, if 
one were to have Mj — M m for j ^ to, then 9j = 9 m , and 
we showed above that identical bare momenta results in 
the Bethe ansatz wavefunction vanishing. Summing l|47|) 
over j and using the fact that L(9) = —L(—9) means 
that the total momentum is 



(48) 



The 9j in the ground state are all real, and lie in the 
region \9j\ < A, where A depends on the chemical poten- 
tial h. There is one particle in the ground state for each 
solution of (|47|l with \9j\ < A; removing any particle in 
this region increases the energy. 

To make any more progress, one needs to go to the 
thermodynamic limit, where there are many particles, 
and the rapidities are closely spaced together. We then 
define the density of particles p(9) so that Np(9)d9 gives 
the number of particles in the ground state with rapidities 
between 9 and 9 + d9. In the ground state, there is a 
particle associated with each integer Mj, as long as the 
rapidity is below A. For \9\ < A, the equations l)47|) 
therefore become 



2np(9) = Rk'{9) 



d9' $(9 - 9')p{9') 



(49) 



-A 



where $(#) = L'(9). The integral equation (gHJl deter- 
mines the ground-state density p(9) for a particular A. 
The maximum rapidity A depends on h p , and is deter- 
mined by minimizing the energy (|44|l . subject to p(9) 
obeying l|49() . In the thermodynamic limit, the ground- 
state energy is 



E p = N 



d9 p(9)e (9) 



(50) 



where the bare energy of a kink is 



e o (0) = —hp - A 



sin 2 (7) 



cosh(2#) — cos(7) 
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The minimum value of the last term in eo is — (1+cos 7) = 
A — 1, so if h p < — 1, it becomes energetically favorable 
to have an empty ground state. Converting back to the 
original variables, this is the Pokrovsky-Talapov transi- 
tion occuring at a = —13/9. 

Even though the density p(9) and the largest rapidity 
A completely characterize the system at zero tempera- 
ture, extracting the Luttinger parameter requires further 
work. We first compute the velocity v of the excitations. 
We define E\ w \ c {9h) an d ^hoie(^h) as the energy and mo- 
mentum change in the system resulting from removing a 
particle of rapidity 9h from the ground state. The veloc- 
ity is then 



dEk 



hole 



<9fchoi( 



dEh \ c /d8h 



9 h =A 



dk ho i c /d9h 



Say we remove the particle associated with integer Mh 
and rapidity 9h- Because the particles are coupled, the 
momenta of all of them changes when a particle is re- 
moved. However, the integers Mj do not change, so from 
gBJ, we see that k hoic = -{2M h + l)n/N. Using gZJ, we 
can rewrite this in the thermodynamic limit in terms of 
the rapidities. After a few manipulations, we have 



dk 



hole 



d6 h 



-2irp{9 h ) 



(51) 



just as in the case without hard cores. 
The energy of a hole of rapidity 9h is 



hole 



= -eo( 



where 56 j is the change in rapidity of particle j when the 
hole is created. Because the integers Mj do not change, 
(|T7j) requires that 

= RNk'iO^SQj + \k{0j) - Afchoie - L(0j - 9 h ) 

rn^h 

Taking the thermodynamic limit and using (|49(l simplifies 
this to an integral equation: 

2nNp(9)86{9-9 h ) = L(9 - 9 h ) - \k(9) + \k holc (0 h ) 

-N [ d9'<S>(9-9')p(9 , )59(9 , ;9 h ). 
J-A 

Note that 59 is a function of both 9 and 9^. Let us define 
the function t(9) as the solution of the integral equation 

i-A 



2itt(9) = -$(0 - A) + 



d6'$( 



')r{9'). (52) 



because fchoie(#/i) is odd in Oh- Because e' (9) is odd in 9, 
we use this with expression for hole energy to get 



dEk 



hole 



-4(A) 



-A 



d9 c' {9)t(9) (53) 



As with fchoio, this is the same as in the case without 
hard cores. Combining this with l|51|) gives the velocity 
of the excitations. We can check this by noting that as 
a — » —13/9, there are no p bosons in the ground state. 
This means that A — » 0, so it follows from (|53|) . H51fl and 
(l4*9l that i! — > e' (0)/k'(0) = 0. The quasiparticles here 
have dispersion i^hoio on fchoio' as expected at a z = 2 
Pokrovsky-Talapov transition. 

The computation of the Luttinger parameter l|46(l is 
affected by the hard-exclusion effects. In particular, be- 
cause R = (1 + N p /N)/3 in l|49|) depends on the number 
of particles N p , dR/dA 7^ 0. First we compute how the 
number of particles changes when A is varied. It is con- 
venient to define the rescaled density p — p/R, so that 
for fixed A, p is independent of R. Then if we define x(#) 
as the solution of the integral equation 

,A 

2w X (9) = L(9-A)+L(9+A)+ / d9'$(9-9') x (9') (54) 

J-A 

it is straightforward to show that 



dp(9) 
OA 



d X (9) 
1-X(A) 39 ■ 



Because 



we have 



N 



Np/N' 



NL 



dNj 
~d~A 



" - NR* 2?(A) 



1 - X(A) 



To compute the Luttinger parameter, it is convenient 
to also define a rescaled energy E — E p /R. A little 
algebra then yields 



1 d E'R? 



d 2 E v 



ON* RN!p OA \ N p 



where E' = dE/dA. Computing E' is then similar to the 
R = 1 case of [2fjj, yielding 



We then have 
d 



N—[p{9)89(9-9 h )- 
= t{6)-t{- 



p(-9)59(-9;9 h )} 



Since none of the quantities inside the derivative depends 
on R, we can use the calculation of [2(j to get 

d*E p 1 (l-x(A)) 2 ( f A \ 

-QNJ ~ & 2p(A) [ e ° {0) ~ J_ A ^ e)T{e) ) 
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where t{&) is defined above via l|52|) . Combining this 
with our computation of the velocity above yields finally 
the Luttinger parameter 



K„ 



1 + N p /N 
1 - X(A) 



(55) 



The effect of the hard-core constraints ends up being the 
(1 + Np/N) 2 in the numerator. 

We have seen that when the coupling a — —13/9, there 
are no p bosons in the ground state. The equation i|55|) 
gives the correct K p = 1, because x(0) = and N p = 
here. As a is increased, p bosons with \6\ < A fill the 
ground state. For some value of a, A — * oo. We can find 
out how many bosons there are here, because in this limit 
we can solve the equation (|49|) by Fourier transformation. 
Actually, to just get N p (A — > oo), we do not need to go 
to the trouble: we know from the XXZ model that 



lim 



d6p{d) 



-A 



Thus in our case, we have N P (A oo) = N/5. The 
equation l|55|l for the Luttinger parameter is therefore 
valid only for values of chemical potential such that only 
solutions of the Bethe equations with all bare momenta 
real are present in the ground state. In other words, 
(|55|1 is valid only for a such that N p < N/5. To find the 
Luttinger parameter for the remaining range of tr, we can 
repeat the above analysis for the t bosons. This yields 



K t = 



1 - N t /N 
1 - x(A t ) 



(56) 



where A t is the maximum rapidity of a t boson in the 
ground state. By the same argument, this expression is 
valid for N t < N/5. 

Because 2N p + 3N t = N, the expressions for the Lut- 
tinger parameter are both valid at exactly one point: 



where A 



A, 



oo, 



N, 



N t = 1/5. This pro- 
vides a good check on our result: at A = A t = oo, 
N p = Nt = 1/5, comparing the expressions (|55|l and 
(|56[) yields K t — 4K p /9, as noted in general in (|35[1 . In 
fact, we can find the exact value of K p at this point, by 
either a Weiner-Hopf analysis |2l| , or by appealing to the 
XXZ model, where the answer is known. Either way, one 
finds that 

lim (1 - x(A)) 2 = 2cos- 1 (A)/tt 

A-+oo 

For our case of A = 1/3, we have therefore K p = 1.83754, 
so vortices are relevant here. Another check on our re- 
sults is to note that the computation of the velocity v 
must give the same answer at this point whether com- 
puting with the t bosons or the p bosons. Indeed, we 
see from (15 II) and (|53|l that both are independent of R 
and A, only depending on A; since A is the same for both 
pictures, we do get the same answer. 

The incommensurate regime is characterized by having 
Kp < 9/8. We cannot analytically find the value of a 



Luttinger parameter 



P 1.6 




vortices irrelevant 



FIG. 9: Bethe ansatz results for K v as a function of a. Recall 
that K t = 4K p /9. 



where K p — 9/8, but it is straightforward to solve the 
above integral equations numerically. First, one must 
determine the value of A for a given a by solving l|49|) 
for p(9) and A while demanding that the energy l|50|) 
be minimized. Integrating p(9) yields the ground-state 
particle density N p /N as well. Knowing A then allows 
one to solve (I54|l for x(#) and hence x(A). Plugging this 
into H55fl then yields the Luttinger parameter. We find 
that K p = 9/8 when a = 1.422.... This leads to our final 
result for the phase diagram: the incommensurate phase 
exists only for 

1.422 < a < 1.444 . 

The incommensurate region is also very narrow near the 
Potts critical point M3, so it is likely that it occupies 
a very small region of parameter space. We plot K p as 
obtained above in Fig|5] 



VI. CONCLUSION 

The primary purpose of this paper was to provide a 
thorough study of the very simple one-dimensional boson 
model in ©■ This model was originally motivated by 
the studies of 'tilted' Mott insulators of atoms in optical 
lattices. Despite its simplicity, our model has a rich phase 
diagram, shown in FigfJ with a many competing ordered 
phases and possibilities for exact solutions. 

Our solutions shed additional light on the nature of the 
phase transition in the chiral clock model. Over a certain 
regime of parameters, we provide definitive evidence for 
a gapless, floating incommensurate phase adjacent to the 
gapped Z3 ordered phase with a period-3 density wave. 
However, the size of this incommensurate phase was re- 
markably small, making it unlikely that it will ever be 
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observed in experimental or numerical studies. We also 
obtained specific realizations of two of the exactly solu- 
ble multi-critical points identified by Andrews et al. [13 , 
which are now known to be part of the series of minimal 
conformal models 0. It is likely that our analysis could 
be generalized to this series; it would be interesting in 
particular to see if the incommensurate phases existed in 
general. 
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APPENDIX: CONNECTION TO BAXTER'S 
HARD SQUARE MODEL 

In this appendix we show that the quantum Hamilto- 
nian Ji in © generates precisely the transfer matrix of 
the classical, two-dimensional hard-square model intro- 
duced by Baxter Q. 

It is useful to note that the hard-square model, as de- 
fined in the Introduction, can also be interpreted as a cer- 
tain 'interface' or 'height' model. In particular, it is iden- 
tical to the so-called A4 restricted-solid-on-solid (RSOS) 
[23| . To see the connection with the RSOS model, we 
need to associate each tile configuration with a corre- 
sponding set of heights. On one sublattice of the square 
lattice, associate each site with a tile present with height 
1, and a vacancy with height 3. On the other sublat- 
tice, associate each tile with height 4, and a vacancy 
with height 2. Then the hard-square restrictions are sim- 
ply encoded in a set of constraints on heights on nearest 
neighbor sites: height 1 can only be next to 2, 2 can be 
next to 3 or 1, 3 can only be next to 2 or 4, and 4 can only 
be next to 3. These height restrictions are conveniently 
summarized in the A 4 Dynkin diagram, and hence the 
terminology. Like the hard-square case, the general A n 
RSOS models are solvable along two lines, with one crit- 
ical point on each. These critical points are multicriti- 
cal points with known conformal-field-theory descriptions 

Here, we continue to use the terminology of hard 
squares. We consider the diagonal-to-diagonal transfer 
matrix of the hard square model. In the orientation of 
Fig. ^2 this transfer matrix acts on the space of states 
of a zig-zag line of lattice sites as shown in Fig. 1101 

There are three interaction parameters: L and M are 
the interaction strength of adjacent tiles, while z is the 
fugacity for each tile. The grand canonical partition 
function can then be written as the sum of products of 




FIG. 10: Zig-zag line defining the space of states upon which 
TC acts. 

Boltzmann weights. The Boltzmann weight of a plaque- 
tte is shown in Fig. 1111 where the heights a, b, c and d 

c 

I) d _ ^(a + b+c+d) / i^Lac+Mbd^ — a+b—c+d 

a 

FIG. 11: Boltzmann weight of a plaquette of the square lat- 
tice. The tiles are centered on the vertices of this plaquette, 
and a, b, c, d = 0, 1 indicates absence/presence of a tile. 

are or 1; 1 denotes a tile, so the restriction is that 
ab = be = cd = da = 0. The extra parameter t cancels 
out of the partition function and so is arbitrary. The pure 
hard square model is obtained by setting L = M = 0, 
but we will need the full range of values of these cou- 
plings to explore the parameter space of Ji. The hard- 
hexagon model (on a triangular lattice) is obtained by 
setting L — and M — * — 00. 

We will now take the so-called r-continuum limit of 
this transfer matrix, so that we may obtain the corre- 
sponding Hamiltonian which evolves the states in a con- 
tinuous imaginary time. To take this Hamiltonian limit, 
we define t = exp(i/4) and z = e~ L (l — £). Then 
the Boltzmann weights for all the allowed possible states 
around each plaquette are shown in Fig. ^| 

When ( = Al = and L — > 00, the vertical trans- 
fer matrix is the identity. The Hamiltonian is found by 
expanding around this limit. The small parameters are 
then C, M and e~ L l 2 . Note that to get back to the orig- 
inal Hilbert space, we need to break the transfer matrix 
into two pieces. The first evolves the lower of the two 
rows of sites in the above zig-pattern, and then the sec- 
ond evolves the upper row. The full transfer matrix is 
then 

TiT 2 « (1 - CWi)(l - (H 2 ) « (1 - (H), (A.l) 

where TL = Hi + TL2 ■ Let d| be a boson creation opera- 
tor which creates a tile on site j. We identify these tiles 
with the bosons of 0. It is easy to see that the con- 
straints of the hard square model correspond precisely to 
the constraints J2J. Moreover, here we obtain the same 
Hamiltonian in (up to an arbitrary overall energy 
scale) with the couplings 

-=e L ' 2 C ■ ^=-Me L / 2 (A.2) 
w w 

Note that M > means attractive next-nearest-neighbor 
interactions. 
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(1 - C) 1/2 e M 




1/2 




o = (l-c) 



i = (i-c) 1/4 




l/4„-£/2 



= (l-Cf^e 

1 


= 1 



FIG. 12: Boltzmann weights for the possible states of a pla- 
quette of the hard square model 



Baxter Q showed that the hard-square model is inte- 
grable for 



(l-e- L ){l-e~ M ) 



(A.3) 



In the Hamiltonian limit, this condition is equivalent to 
Q . Baxter also showed that there are two critical points 
within the parameter space of i|A.3|) . Defining 



I = z- 1/2 (l-ze^ +M ), 



L+M\ 



the critical points are at 



I = ± 



1 + V5 



-5/2 



which leads to QfilSjl . The critical point at © is in the 
universality class of the tricritical Ising model, while that 
at |JSJ is of the 3-state Potts type jj Q . The physical 
interpretation of these critical points is discussed in Sec- 
tion ITU 
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